Mesoscopic phonon transmission through a nanowire-bulk contact 
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We calculate the frequency-dependent mesoscopic acoustic phonon transmission probability 
through the abrupt junction between a semi- infinite, one-dimensional cylindrical quantum wire and 
a three-dimensional bulk insulator, using a perturbative technique that is valid at low frequency. 
The system is described using elasticity theory, and traction-free boundary conditions are applied 
to all free surfaces. In the low-frequency limit the transmission probability vanishes as lj^, the trans- 
port being dominated by the longitudinal channel, which produces a monopole source of elastic 
radiation at the surface of the bulk solid. The thermal conductance between an equilibrated wire 
nonadiabatically coupled to a bulk insulator should therefore vanish with temperature as T^. 

PACS numbers: 63.22.-|-m, 85.85.-|-j 
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I. INTRODUCTION 

Electronic transport in a variety mesoscopic systems 
has been successfully described by the theory of Landauer 
and Biittikeri ^i^i?i'^i^i^i7 In this scattering approach, each 
fully propagating channel in a wire contributes 27re^/?i 
to the electrical conductance. Recently there have been 
experimental efforts to study phonons in the mesoscopic 
regime^ and beautiful experiments by Schwab et al^ and 
by Yung et aliS, have observed the low-temperature uni- 
versal thermal conductance quantum of Trk^T/Gh per vi- 
brational channel. In these experiments, the nanowircs 
were adiabatically (on the scale of the thermal wave- 
length) connected to thermalizcd phonon reservoirs^ 
and the observation of a thermal conductance varying 
linearly with temperature implies that very little phonon 
reflection occurred at the nanowire-bulk junctions. 

In the presence of an abrupt, nonadiabiatic cou- 
pling between the nanowire and bulk reservoirs, phonons 
will scatter at the junctions and suppress the thermal 
conductanceAiSiiiii Cross and Lifshitai^ have calcu- 
lated the frequency-dependent acoustic phonon transmis- 
sion probability T between a semi-infinite quantum wire 
of rectangular cross-section and a thin plate with the 
same thickness as the wire, and find that T cx w^/^ in the 
low-frequency limit. In Ref. |l4l | a short nanowire, mod- 
eled as a harmonic spring, abruptly connected at both 
ends to three-dimensional bulk insulators was found to 
have T (X uj^. These result suggest that such nanowircs 
will eventually become thermal insulators at low temper- 
atures, and the universal thermal conductance quantum 
will not be observed. However, such a cross-over to insu- 
lating behavior has not yet been observed experimentally. 

Understanding the scattering caused by nonadiabatic 
nanowire-bulk contacts will be essential for pushing 
phonon physics into the mesoscopic regime, as well as 
for the design and operation of thermal nanosensors such 
as calorimeters and bolometers. In this paper, we extend 
the previous work by calculating the mesoscopic acous- 
tic phonon transmission probability through the abrupt 
junction between a semi-infinite, one-dimensional cylin- 
drical quantum wire and a three-dimensional bulk insula- 



tor, using the low-frequency perturbative technique em- 
ployed by Cross and Lifshitz:^"- xhe nanowire and bulk 
insulators arc described using isotropic elasticity theory, 
and traction-free boundary conditions are applied to all 
free surfaces. In the low- frequency limit the transmis- 
sion probability is found to vanish as o;^, with the low- 
frequency transport being dominated by the longitudinal 
channel. 

In the Sec. we review the calculation of the long- 
wavelength vibrational modes of an infinitely long cylin- 
drical elastic rod. The long wavelength limit of interest 
here is defined as /c6 <C 1, where b is the radius of the 
cylinder and k is the wavenumber. In Sec. IIIII we show 
that in the long- wavelength limit, the bulk solid produces 
a hard-wall boundary condition on the nanowire. In 
Sec. IIVI wc calculate the displacement field in the three- 
dimensional bulk solid given an applied traction to its 
surface, using what is essentially an elastic Green's func- 
tion mcthod^^^ The frequency-dependent transmission 
probabilities for each of the four gapless modes are calcu- 
lated in Sec.0 and our conclusions arc given in Sec. lVIl 



II. 



VIBRATIONAL MODES OF CYLINDRICAL 
WIRE 



In this section, we will briefly review the elastic waves 
for an infinitely long cylindrical waveguide. In the long 
wavelength limit there are four branches, which include 
one torsional branch, one longitudinal branch, and two 
flexural branches 4^ We assign a numerical subscript to 
represent each branch mode, with "1" denoting the tor- 
sional branch, "2" denoting the longitudinal branch, and 
"3" and "4" denoting the flexural branches. These four 
branches are orthogonal. Also, cylindrical coordinates 
are used below. 

We assume an isotropic elastic continuum with trans- 
verse and the longitudinal sound speeds 



and 



'A + 2/^ 



(1) 



where p is the mass density, and A and fi are the Lame 
constants. 



A. Branch 1: torsional 



The displacement field is given by 

ui(r,t) = r e'C^^-^i^'ee, 



with dispersion relation 



oji = Ctk. 



(2) 



(3) 



The stress tensor elements acting on a cross-section of 
the rod are 



fywz — M 



\ dz dr 
f duie 1 dui 
\ dz r 86 



A(V-Ui) + 2/i 



dz 



= 0. 



(4) 



B. Branch 2: longitudinal 

The displacement field is 

U2(r,t) = [/,(r)e, + /,(r)e,] e^C^--^*), 

where 

jr(r) = -A2aJi{ar) + B2ikJi{f3r), 
f,{r) = A2ikJo{ar) - B2/3Mpr). 



(5) 



(6) 



J„(r) is the Bcssel function of nth order, a and (3 are 
constants determined by the boundary conditions. Fur- 
thermore, 



and 



A2 ^ f3^~k^ JiiPb) 
B2 2iak Ji (ab) 

The relevant stresses are 



(7) 



(8) 



(^2rz = M 
0'2ez — 0, 



^22 



A 



ikfr{r) + 

dr 

dfrjr) ^ frjr) 
dr r 



i{kz—uj2t) 



+ifc(A-f 2M)/.(r) 



^i{kz—uJ2t) 



(9) 



In the long wavelength fc6 <C 1 limit, lu2 = cgk, 



a = ik^l-^, and p = k^^~l, (10) 
where cq is related to Young's modulus of elasticity E by 



To leading order 

U2 (r,i) =e, eH^-— 
and the stresses are 



(11) 



(12) 



0'2rz 


= 0, 


0'2ez 


= 0, 










0'22z 


= lii — -k 















where p = ct /ci is the ratio of the transverse to longitu- 
dinal sound speed. 



C. Branch 3: i-polarized flexural 

The displacement field is 



U3(r,t) 



gr (r) cos e,. + gg (r) sin 9 eg 



-f- gz{r) cos Ogz 



J{kz-uj3t) 



(14) 



where 

9r{r) 
9e {r) 



dMar) , dJi(^ Ji(/3r) 

- B3 ; h O3 , 



dr 



dr 



Ji(ar) Ji(/?r) dJi(/?r) 
i^a O3- 



dr 



g^[r) = ikJi{ar) - iB3 — Ji{Pr). 



(15) 



B3 and C3 are constants. The stresses are given by 



ikgr{r) 
ikgg{r) 



dgzjr) 

dr 
9z{r) 



sine e'^^'-'^^^\ 



A 



dgr{r) gr{r) + geir) 



dr 

i {X + 2fi) k g^{r) 



cos 



^gi(fcz-a.3t)^ (16) 
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In the kb << 1 limit, 



^3 = yofc , 



a = ik\ / 1 



/ cpbk 
\ 2ci 



and 



(17) 



(18) 



(19) 



To leading order, 

U3(r,i) = (cos ee,. -sin 6*66 + ifca;e^)e'('''^-'^^*' 
= (e^-i/cxe^)e'('=^-"^*), 

and the stresses are 

2 

]_ 
T 



(20) 



Q '^0 \ 2 



2 



The rod bends in the xz plane 



(21) 



D. Branch 4: y-polarized flexural 

An independent flexural mode can be found by letting 
the rod bend in the y direction. In the long-wavelength 
limit, 



presents a stiff boundary to the nanowire, so to zeroth 
order one calculates the vibrational modes of the wire 
assuming a zero-displacement boundary condition at the 
contact. The stress distributions associated with these vi- 
brational modes is then calculated in the junction region. 
These zerotli-order vibrational modes, however, do not 
carry elastic energy because of the hard-wall boundary 
condition associated with the infinitely stiff bulk solid. 
Now one relaxes the hard-wall boundary condition, re- 
placing it with the condition that elastic waves in the 
bulk are purely radiative, having outward moving com- 
ponents but no inward ones. The elastic energy radi- 
ated by the nanowire into the semi-infinite bulk solid is 
then computed using the actual elastic parameters of the 
bulk, and the ratio of incident to transmitted energy de- 
termines the transmission probability. 

Thus, as explained, we will calculate the elastic stress 
on the surface of the three-dimensional bulk insulator, 
produced by the vibrating nanowire, by assuming that 
the bulk imposes a zero-displacement boundary condition 
on the nanowire. This is physically reasonable, and can 
be further justified by considering the bulk to be a thick 
wire with a radius B much larger than the nanowire ra- 
dius Assuming a sound wavelength larger than both 
b and B, the conservation of linear and angular momen- 
tum lead the zero-displacement boundary condition in 
the hmit i? > 6. 

We consider a semi-infinite cylindrical elastic nanowire 
of radius &, lying along the z axis and connected at 2; = 
to a thicker cylinder of radius B. An incident elastic wave 
Wi ~ \ii (r, 9) e'C^^""^*) propagates from the nanowire to 
thick cylinder, k is smaller than both h and B so both 
cylinders are still one-dimensional, and i = 1,2,3,4 de- 
notes the branch. The scattering causes a refiected wave 
for z < and a transmitted wave for z > 0. We can write 
the displacement field as 

u, (r, 9) e^C^^-"*) + R,ju* (r, 61) e-*('=^+'^*) for z < 



U4 (r,i) = (e,y-ifcyez)e»('=^-"^*\ 
and the stresses on the surface normal to z are 



(22) 



CTiz 



T 



3 



(23) 



r„-Uj(r, 9) e 



(24) 

Here Rij and Ty are the refiection and transmission co- 
efficients, which are matrices in the channel indices. 

The continuity of the displacement field, combined 
with the orthogonality of the vibrational eigenmodes, 
leads to 



Rij — Tij . 



(25) 



In the Appendix we show that in the B ^ b limit, con- 
servation of linear and angular momentum leads to 



R, 



—Sji and Til 



0, 



(26) 



III. BOUNDARY CONDITIONS AT THE 
NANOWIRE-BULK INTERFACE 

The essence of the pcrturbative method we use is as 
follows: In an abrupt junction geometry, the bulk solid 



which means that elastic waves are flipped upon reflec- 
tion, and no interbranch scattering occurs. This result is 
analogous to that obtained by Cross and Lifshitz in their 
thin-plate geometry^ 

Linear combinations of the vibrational eigenfunctions 
of Sec. |n] can be used to satisfy the boundary conditions 
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of Eq. (|26() . These displacement fields produce the fol- 
lowing stress distributions at the z = interface: For the 
torsional mode we obtain 



2i^fcHe~*"* for r < b, 



for r > b. 



(27) 



For the longitudinal mode, 



2i/zf^fce-"^* for r < b, 



(28) 



for r > b. 

And for the x-polarized flexural mode, we find 

(Jxz (^yz ^ (^zz ^0 for r > b, (29) 



(l + ^) {b^ - x^) 



3 - ^ ) 2/2 



for r < b. 



(30) 



i/i ( 1 - 1^ ) k^xye 







The stress distribution from the y-polarized flexural 
mode has the same form as ()29|l and (|3U|I after exchang- 
ing X ^ y. 



IV. 3D ELASTIC RESPONSE FUNCTION 

Next we calculate the displacement field in the three- 
dimensional solid given the applied stress of Sec. IIIII to 
its surface at z = 0. For r < b, this is the stress distri- 
bution produced by the nanowire, and for r > & it is the 
stress imposed by the traction-free boundary condition. 
The method we use here, which is essentially a Green's 
function method, is well known in elasticity theorypiii^ 

To find the displacement field in the bulk solid given 
the boundary conditions described above, a scalar poten- 
tial (/) and a vector potential H are introduced according 
to 

u = V(f> + V xU. (31) 
The wave equations for the potential fields are 



^-cM^ = 0, (^-C.2V2)H = 0. (32) 



They can be written as 
1 C'^ 

cl) (x, y,^) = ^ dhdhf (fci, fe) e-'^^^^+'^^y^e''^^'' 



(33) 



and 

H {x, y, z) 
where 

and 



1 

2^ 



c?fcic?fc2h (fci, e 



ki3 



'i{kix~{-k2y) ikt3Z 

(34) 
(35) 



Ct3 



(36) 



Here / and h are the Fourier transforms of the potential 
fields (j) and H at z = 0. 

Now, we can use Eq. (|31|l and choose the transverse 
"gauge" V • H = to express the components of the dis- 
placement vector by the inverse Fourier transform J-~^, 

{x, y, z) = -iT-^ [fci/ (fci, fc2) e*-^^ + g., (k^M) e**^'""] 
Uy (x, y, z) = -iT-^ [fc2/ (fci, fc2) e*''-^' + (fci, fe) e*'='^'] 

fci3/(fci,fc2)e*-^^ 



U2 {x,y,z) = iT 



fclffa: (fcl,fc2) + fc2ffj/ (fcl,fc2) ifct: 



(37) 



where 



5a;(fcl,fc2) = fc2^2 (fci, fc2) + fctS^y (fcl,fc2) (38) 



9y (fcl,fc2) 



-kih^{ki,k2) ~ kt3hx{ki,k2) ■ (39) 



The stress at the boundary z = can also be expressed 
in terms of the inverse Fourier transforms, as 



dux du. 



dx 



dz y 
^iT-^ 2kiki3f{k,,k2) 

{hz^ - fcl^) gx (fci, fc2) - fclfc25y (fcl,fc2) 



«t3 



duy du^ 
dz 

-1 



(fcta^ 



9y / z=o 

2fc"2fc;3/(fci,fc2) 
- k2'^) Qy (fci, fc2) - kik2gx (fci, fc2) 



«t3 



A (V • w) + 2^i^ 
dz 



= iiT-^ [(fci^ 



fc2 



z=0 



4*3 



/(fcl,fc2) 



+ Sfciffa; (fcl,fc2) + 2fc2.gy (fci, fc2)] 



(40) 



By giving the boundary values of Gxz, ^^yz, and Gzz, 
we can find /, gx^ and gy from the equations above. If at 
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least one of these three stresses is nonzero, we obtain 



1 J 




+ - kt3^ 


)^ 




(41) 


Vo (fci, ^2) 


A* 


+ 2kt3hT 






Oyz 


}■ 


(42) 



9x 



Vo iki,k2) 
"771 (fci,A:2) 



2ki3kt3kiT 



kt3 
m (fci,fe) 



iktski" 



■^13 



2kt3kik2 


T 


Oyz 


} 











and 



1 



9v 



r\Q (fci, ^2) 
'r\2 (^2, fci 



kt3 
r\\ (fc2, ki) 



2fc(3fct3fc2^ 

2kt3k\k2 



k 



t3 



2fct3fe^ 






} 











(43) 
(44) 
(45) 

(46) 
(47) 
(48) 



Here 



Vl{ki,k2) = k2^ + ki^k2^+ki^kt3^ 

+ 2kt3k2^2ki3 - kt3) + kt3^ (49) 

m{kuk2) = fcifca [fci' + fc2' + fct3(4A:/3-3fct3)](50) 

ki 

?7o(fci,fe) = miki,k2) + —i']2{ki,k2) . (51) 

k2 

Therefore, from Eq. ()37|l . we can find the displacement 
vector at z = in terms of the boundary stresses. 



= 



kt3'no (fci, ^2) 



+771 {ki,k2)T 
^^1 



fct3% (^1,^2) 



-?7i {k2,ki)T 



kikt3ri3 {ki,k2)T 

- m {ki,k2)T 
k2kt3m {ki, k2)T 

- m {k2,ki)T 





}] 







C^xz 


}] 







'no{ki,k2) [ct 



—ka^ 



-m (fci, ^2) kiT 



k2T 



V. ENERGY TRANSMISSION FROM 
NANOWIRE TO BULK 

Now we are ready to calculate the transmission prob- 
ability, defined as the ratio of transmitted to incident 
elastic energy flux, for each of the four gapless branches. 
The energy current / can be expressed asii 



W"^' ' dt 



I ^ { I dxdy I -gf^xz + -KT^^yz + -^<^z 



z=0 



where (...) represents a time average and J^dxdy is the 
surface integral over the z = cross-section of the wire. 
In conventional complex notation for waves, this becomes 



/ = -| Im ^ dxdy {u^ct'*^ + UyCr'*^ + u^a^^ . (55) 

We will calculate the energy current for the different 
branches separately. 



A. Torsional branch 

The stress distribution (|27|1 in rectangular coordinates 

is, for x"^ + y"^ < b, 



(Jzz = 

o-xz = -2i^fcye"*'^* 
ayz = 2i^kxe~^'^^ , 



(56) 



and zero for y/x'^ +y'^ > b. Using Eq. we obtain 



-^^ Im^ dxdy {u^a*^ + Uya*y^) 

Re / '^^^/^^ ^ / d.dye-'(^-+^^^' 

47r J fct37?0 (fcl, fc2) 7s 



(57) 



X {ki,k2)T 
+ (^*yz [rii {k2,ki)T 
^ 2/ifc^w Re 



0-:r 



yz 



dkidk2 



- m {ki,k2)T 
m {k2,ki)T 



kt3-qa (fci, fe) 
X {T[y]* {ini{ki,k2)J'[y]+m{ki,k2)T[x\) 

+ F[x\* {in{k2,ki)F[x\+r^2{k2,ki)T[y])] . (58) 

Expanding the Fourier transforms J- [x\ and T [y] for 
small fc6, and keeping only the leading terms, 



(52) 



with 



(59) 



7?3 {k^,k2) = fci' + k2^ + kt3 {2ki3 ~ kt3) . (53) 



2tt _ ^ 

1_ I y^-^ik.x+k..y)a^ay (60) 
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we have 



Flexural branches 



32 

X 



Re 



dkidk2 



fci 



771 (fci,fc2) + — 7?2 
^2 

Hb^k'^uj _^ /■ „ „ ^2^ + fci^ 



32 



Re / dkidk2- 



kt3 



24^ 



(61) 

(62) 
(63) 



By normalizing / with the energy current j^b^kuj of 
the incident torsional wave, we obtain the transmission 
probability 



T = -b^k\ 
6 



B. Longitudinal branch 
Using Eq. if^ . 

2ti 



(64) 



Jo 



Im / rdr dO {u^al^) . (65) 



Then to leading order in kb, 



b i'27T 

rdr / d9(7:,^ 

"'0 



Re / dkidk2 



uj'^kis 



ct^m (^1,^2) 



= ^{%l,^ku] Ro I dk,dk2- "'^'^ 



2ct V ct 



'ctT]o iki,k2)' 



(66) 



Normalizing by the power j^^b kuj carried by the in- 
cident wave leads to 

^(fc) = ^&'fcc^Re / dkidk2 '^.t" , , (67) 

TTCt'^ J Ctr]o{ki,k2) 

= Ub^k^, (68) 



where 



t; = ° „ Re / dkidk2- ''^ 



■ct77o {ki,k2)' 



(69) 



Assuming the materials to be made of Si, we have p 
0.694 and 



r = 1.916^fc^ 

0.923 f|^y 



(70) 
(71) 



By using flexural stress distribution H30() . we find that 
to leading order in kb the energy current is 



— Re / dxdyT 
^ J s 



fiLok^ 
32^ 
2 
2" 



axzdxdy 



m (fci,fc2) 
fct3?7o (fci,fe 

Re y dkidk2 
2 



^yi (fci, ^2 



Co „ » 2 

, - 3 2/ 



fct3??0 {kl,k2) 

dxdy (^l + - x^) 

Vi (fci,fe) 



Re / dkidk2 



ktsm [ki, ^2) 



^ /^"^O ,4,3\ T, /j, J, Vi{ki,k2) , . 

= ^ — ^ w Re / dkidk2- t.A -(72) 

32 Vct"^ / J fct3?7o (fci, k2) 

Normalizing by the energy current of the incident wave 
/in = Im J^dxdy {a* ^ - ikxal^) ^^^Q (73) 

(74) 



^^^0 ,4,3 
4 ct^ 



where we have used Eqs. (|20|l and H21|) . the transmission 
probability is found to be 



tfb^k^ 



b^k^uj Re / dkidk2 



CtTIi (fcl,fc2) 

'i^ktsVo {ki,k2) 



where 



Co 



Re / dkidk2 



Cf7?i (fci,fc2) 

wfcta^ (fci,fe) 



(75) 
(76) 



(77) 



Using p ~ 0.694 for Si, the transmission probability 
becomes 



T = 0.268 



0.609 



5/2 



(78) 
(79) 



Because of the cylindrical symmetry, the y-polarized flex- 
ural branch has the same transmission probability as the 
x-polarized branch. 



VI. CONCLUSIONS 

On the left side of Table d we summarize the trans- 
mission probability results calculated above, as well as 
the low-frequency dispersion relations of the four gapless 
modes. For comparison with the results of Cross and 
Lifshitai^ for a rectangular wire connected to a thin plate, 
we reproduce their results on the right side of this Table. 
In each case there are four gapless acoustic modes: one 
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torsional, one longitudinal (or compressional) and two 
flexural bending modes. Also, the form of the dispersion 
relations are the same for both wires. For all branches 
the transmission probability to a three-dimensional bulk 
solid has a higher-order frequency dependence. This is 
at least partially a consequence of the higher vibrational 
density of states in the three-dimensional system as com- 
pared to a plate: For the longitudinal and x-polarized 
flexural branches, T is one order higher in w, consistent 
with the density of states enhancement 

The phonon transmission probabilities can be used to 
calculate the mesoscopic thermal conductance between 
an equilibrated wire and bulk. According to the thermal 
Landauer formula^^'ii'^ a total transmission probabil- 
ity T{uj) varying at low frequency as to'^ will lead to a 
low-temperature thermal conductance varying with tem- 



perature as Gth oc T'^+i. In our case, T(a;) is a sum of the 
T for each channel. The thermal conductance between 
an equilibrated cylindrical wire nonadiabatically coupled 
to a bulk solid should therefore vanish with temperature 
as T^. 

Finally we comment on the applicability of our results 
to nanoscale phonon experiments, which do not consider 
infinitely long wires and perfectly sharp corners. For our 
theory to be valid, the wire must be longer than the sound 
wavelength, and the characteristic radii of curvature at 
the junction must be much smaller than this wavelength. 
Therefore, because of the first condition, our results will 
become invalid in the extreme low-temperature limit, and 
the conductance will crossover from our predicted 
scaling to some other behavior. 



TABLE I: (left) Dispersion relations uj (k) of the low-frequency vibrational modes in a cylindrical nanowire, and transmission 
probabilities T through the junction with a three-dimensional bulk insulator, as a function of both k and uj. ti and tf are 
constants defined in Eqs. 1691 and 17711 . In the low- frequency limit the total transmission probability vanishes as uj^, the 
transport being dominated by the longitudinal channel, (right) Same quantities for a rectangular wire connected to a thin 
plate, reproduced from Ref. H^ . Here Ii and I2 are Poisson-ratio-dependent numbers. 



cylindrical nanowire (radius b) 
— > semi-infinite space (3D solid) 


rectangular nanowire (width 6, thickness d) 
— > thin plate (2D plate of thickness d) 


branch 


uj{k) 


T{k) 




branch 


tj(fc) 


T{k) 


T{co) 


torsional 


Ctk 




I (if) 


torsional 


fctk 


hbk 


^ 2d ct 


longitudinal 


Cok 


ti ibkf 


'.(ft)' 


compressional 


Cok 


4bk 


CO 


flexural (x-direction bending) 


Icobk" 


tf ibkf 




in-plane bending 




i {bkf 




flexural (j/-direction bending) 


same as x direction 


normal-plane bending 




hbk 
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APPENDIX A: MOMENTUM CONSERVATION AND HARD- WALL BOUNDARY CONDITION 



Here we use linear and angular momentum conservation to derive Eq. H26|l in the B ^ b limit. First we equate the 
torques produced by the thin and thick wires on each other. From the rotational stress 

(^ez = A* \ TTTT (Al) 



dz r do 

a,ez{r,9)e''^'''-^'^ + Rija*g,{r,9)e-''^'''+'^*'> for z < 0, 
Ty O-je^(r,6»)e*('=^-'^*) for z > 0, 



(A2) 



the torsional torque t = rdrdOragz is found to be 

f 6^ ((5,16*''^ - Rae-''''') /i/ce"''^* for 2 = 0" 
f B^TaMfce^^'^^-"*) for z ^ 0+. 



Only torsional mode contributes to the torque. 
By equating these torques we have 



Combining this result with Eq. H25|l leads to 



R^l = - 0, i^l 



(A3) 



(A4) 
(A5) 



B*_b^ 



Til = 



2b^ 



Taking the B/b ^ oo limit then leads to Rn — > —1 and Tu 0. 

An analogous result for Ri2 and Ti2 can be derived by considering the total force in z direction, 



' inb^ {5i2e'^'' - R^2e-'^'') f^fike-'"^* for z < 0, 



which gives 

6' i5^2 - R.2) - B^T,2, 

since only longitudinal mode has a nonzero Fz . Then 

R^■2 =T,2=0, 1 + 2 



for z > 0, 



(A6) 



(A7) 
(A8) 

(A9) 
(AlO) 



R22 
T22 



B^-b^ 
B^ + b^ 
2&2 



-1 



0, 



(All) 
(A12) 



for B/b 00. 

By further considering the conservation of momentum in the x and y directions, it is not difficult to derive the 
result quoted in Eq. (P^ . 
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